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ABSTRACT 

A controversy has developed regarding the stellar wind mass loss rates in 0- 
stars. The current consensus is that these winds may be clumped which implies 
that all previously derived mass loss rates using density-squared diagnostics are 
overestimated by a factor of ~ 2. However, arguments based on FUSE observa- 
tions of the P v resonance line doublet suggest that these rates should be smaller 
by another order of magnitude, provided that P V is the dominant phosphorous 
ion among these stars. Although a large mass loss rate reduction would have 
a range of undesirable consequences, it does provide a straightforward expla- 
nation of the unexpected symmetric and un-shifted X-ray emission line profiles 
observed in high energy resolution spectra. But acceptance of such a large reduc- 
tion then leads to a contradiction with an important observed X-ray property: 
the correlation between He-like ion source radii and their equivalent X-ray con- 
tinuum optical depth unity radii. Here we examine the phosphorous ionization 
balance since the P v fractional abundance, q (P v), is fundamental to under- 
standing the magnitude of this mass loss reduction. We find that strong 'XUV" 
emission lines in the He II Lyman continuum can significantly reduce g (P v). 
Furthermore, owing to the unique energy distribution of these XUV lines, there 
is a negligible impact on the S v fractional abundance (a key component in the 
FUSE mass loss argument). We conclude that large reductions in 0-star mass 
loss rates are not required, and the X-ray optical depth unity relation remains 
valid. 



Subject headings: stars: early-type — stars: mass-loss — stars: winds, outflows 
— X-rays: stars 
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1. Introduction 

Over the past several years, questions concerning the validity of what we refer to as the 
"traditional" O-star mass loss rates, M, have arisen owing to the predictions of clumped 
wind models (see recent Potsdam Workshop, Hamann et al. 2008, and references therein). 
Abbott et al. (1981) showed that inhomogeneous winds lead to an enhancement in density- 
squared emission processes such as Ha, infrared and radio free-free, which means that the 
inferred M would be overestimated (see also Puis et al. 2006). Although a clear picture 
of these clumpy structures is still being developed, the so-called "clumping factor" (f c i) is 
believed to be ~ 4 to 5, and the reduction in M scales with \f] cX (see General Discussion in 
Hamann et al. 2008). 

The derived M from diagnostics that are linearly dependent on density (e.g., analyses of 
unsaturated UV resonance line profiles) are expected to be independent of clumping effects 
(Puis et al. 2006) and supposedly should provide more reliable M values. In fact, analyses of 
the P V resonance line doublet (AA 1118, 1128 A) obtained from FUSE observations led to 
the conclusion that traditional M are overestimated by a factor of 10 or more (Massa et al. 
2003, hereafter M03; Fullerton et al. 2006; hereafter F06). Such reductions have far-reaching 
consequences. For example, Hirschi (2008) concluded that the well known evolutionary 
tracks of massive stars could survive M reductions by a factor of 2, but not by a factor of 
10 or more. 

Resolution of this M problem is also of particular importance with regards to the ob- 
served X-ray emission line properties obtained from Chandra and XMM — Newton obser- 
vations. Waldron & Cassinelli (2007) analyzed the Chandra HETGS X-ray line properties 
for a large number of OB stars, and two of their conclusions are directly relevant to this M 
issue. 1) All of the resolved X-ray emission lines are very broad (i.e., HWHM range of 300 to 
1000 km s _1 ), symmetric, and the majority have minimal line-shifts. Although Waldron & 
Cassinelli (2001) were the first to demonstrate that these X-ray line profiles could in fact be 
easily explained by a significant reduction in M, this was inconsistent with the then accepted 
observed M, and they suggested a clumpy or non-symmetric wind structure as a possible 
explanation. 2) The radial locations of the He-like fir (forbidden, intercombination, reso- 
nance) line sources, as derived from their f /i line ratios, typically range from 1.2 to 10 i?*, 
and these distances are well correlated with their respective stellar wind X-ray continuum 
optical depth unity radii which we shall refer to as the "X-ray continuum optical depth unity 
relation" (hereafter abbreviated as XODUR). This relation implies that the traditional M 
(like those of Vink et al. 2000) must be correct (i.e., within a factor of 2 or so). Nevertheless, 
large M reductions have become a widely used explanation of the X-ray line profile symme- 
try problem (e.g., Kramer et al. 2003; Cohen et al. 2006) because the emission from both 
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the near and far side of a star would be observable in an optically thin wind. The subject 
is not yet resolved since Oskinova et al. (2004, 2006) finds that the symmetry can also be 
explained by accounting for the porosity of clumped (or fragmented) winds without the need 
for a reduction in M, but Owocki & Cohen (2006) present arguments against the porosity 
influence. 

The primary goal of this paper is to examine the effects of an excess of hard radiation 
within the He II Lyman continuum on the ionization equilibrium of phosphorus by utilizing 
only outer shell photoionization processes. This is important since the proposed large reduc- 
tion in M is based entirely on the assumed P V fractional ionization abundance. In addition, 
we also consider the sulfur ionization balance because of two arguments used by M03 and 
F06: 1) since phosphorous and sulfur have overlapping ranges in ionization energy, the dom- 
inant stages of sulfur can be used as surrogates for the corresponding ones of phosphorous, 
and; 2) both P v and S v are likely to be the dominant ionization stages throughout the 
O-star spectral range. 

We propose that the primary source of this excess hard radiation is produced by the 
"XUV" spectral energy band which has long been known to have the capability to produce 
anomalously higher wind ionization stages (e.g., Waldron 1984; MacFarlane et al. 1994; 
Pauldrach et al. 1994, 2001). Based on solar physics studies, the XUV lower energy bound 
seems to be rather loosely defined, but the upper limit appears to be fixed at 124 eV (100 A)@ 
where the XUV upper limit represents the start of the X-ray energy band. For our purposes, 
we adopt an XUV radiation energy band defined as 54.4 eV (He II edge) to 124 eV, since 
below the He II edge, the radiation is dominated by photospheric emission. Although XUV 
radiation is un-detectable in O-stars, its effects can be manifested by studying other spectral 
bands. In particular, with regards to the two key arguments used by M03 and F06, we show 
that XUV radiation produces dissimilar changes in the phosphorous and sulfur ionization 
equilibria. 



2. Importance of XUV Radiation 

In this section we use graphical arguments to illustrate the importance of XUV radiation 
with respect to the total radiation cooling curve and the ionization equilibria of phosphorous 
and sulfur. To emphasize the importance of XUV + X-ray radiation, we provide comparisons 
with the case when only X-ray radiation is considered. Although unobservable, we believe an 
observational signature of XUV radiation has already been detected by M03. In their FUSE 



1 Withbroe & Raymond (1984) define a XUV energy range from 25 eV (He I edge) to 124 eV. 
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study of LMC stars, they explicitly state that "...C v is the dominant species for O stars". 
This can only mean that there must be excess XUV emission since the C IV photoionization 
edge lies within the XUV energy range (see Fig. 2). 

We first discuss the expected contribution of XUV radiation to the radiative cooling 
curve (e.g., Cox & Tucker 1969), A(T) = P{T)/(n e n H ) ( erg cm 3 s" 1 ) where P(T) is 
the power per unit volume and n e and uh are the electron and hydrogen number densities 
respectively. Although cooling curves have undergone alterations over the years (e.g., atomic 
data updates and added emission lines), the basic shape of the cooling curve has remained 
intact as shown in Figure 1. This shows A(T) obtained from Raymond & Smith (1977; RS), 
Mewe et al. (1985; MEKAL), and Smith et al. (2001; APED) data. It also shows that the 
XUV contribution to A(T) is clearly important for temperatures between 0.5 and 2.0 MK. 

Since the photoionization edge for P v — > P VI is at 65.03 eV, and that for S v — > 
S VI is nearby at 72.68 eV, it seems plausible that their ionization balances should be quite 
similar. This similarity led to the F06 argument that S v can be used as a surrogate for 
P v. However, as shown in Figure 2, just beyond the P V photoionization edge there exists 
a large collection of intense XUV emission lines (their specific contribution is shown in Fig. 
1) that are located just below the S V ionization energy. Whereas, at energies higher than 
this edge, the XUV emission at this temperature is essentially devoid of emission lines, i.e., 
all lines between the S V edge and 124 eV are at least a factor of 100 times smaller than 
the strongest line just below the S V edge. Consequently, this unique energy distribution 
of XUV lines relative to the energy locations of these photoionization edges indicates that 
the XUV emission should produce different effects on the phosphorous and sulfur ionization 
equilibria. This is illustrated in Figure 3 which shows that the P V photoionization rate is 
10 times larger than the S V rate at the temperature of maximum XUV emission (see Fig. 1). 
Figure 3 includes the C IV rate for comparison, and also shows the expected photoionization 
rates using only the X-ray energy band (> 124 eV). As evident from the displacement of 
these curves, the neglect of XUV radiation leads to underestimates of these rates which are 
substantial at temperatures < 2 MK. This implies that alone, the radiation from the X-ray 
energy band is not expected to have a significant impact on the ionization equilibrium of 
phosphorous, as was recently demonstrated by Krticka & Kubat (2009) in their study of the 
Auger effect on the fractional abundance of P v. 

Our graphical arguments illustrate that the XUV line emission for temperatures between 
0.5 to 2 MK is expected to have a major impact on the ionization structure of phosphorous 
but a relatively minor effect on sulfur. This implies that the sulfur surrogate argument needs 
to be re-examined (see §3). In fact, since the C IV photoionization edge is almost identical to 
the P V edge (see Fig. 2), C IV is a more appropriate surrogate for P v. The main difference 
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is that the C IV photoionization rate is larger than the P V rate (see Fig. 3) due to the 
differences in their cross sections. 



3. Ionization Equilibrium Calculations and Required XUV Radiation 

The stellar wind ionization equilibria of phosphorous and sulfur are calculated in a 
straightforward way to determine the level of XUV+X-ray radiation required to affect the 
fractional ionization abundances, i.e., q (P v) and q (S v). We consider a stellar effective 
temperature (T e ff) range from 27500 to 45000 K which covers the O and early B spectral 
range where P V has been used to study M. The ionization equilibrium is determined 
by adopting an ionization/recombination rate balance approach similar to the one used in 
FASTWIND (Puis et al. 2005). The main differences are: 1) we use the photoionization 
cross sections of Verner & Yakovlev (1995); 2) a wind diffuse field as prescribed by Drew 
(1989), and; 3) we assume a radial power law dependent wind temperature (T w ) which is 
adjusted to produce a phosphorous wind ionization structure similar to that of Puis et al. 
(2008) (i.e., no XUV+X-ray radiation) for all T e // considered. This adjustment leads to a 
T w /T ef f = 1.15(i?*/r) a5 with a minimum value of 0.6T e// . We calculate q (P v) and q (S v) 
for wind conditions at a location where the wind velocity is 1/2 the terminal velocity, v^. 
For each T e ff, a consistent set of stellar parameters (log g, Lb u R*) is found by using the 
fitting-formulae given by Martins et al. (2005). We use the electron scattering Eddington 
factor as described by Lamers (1981) to determine the effective escape velocity, V esc . For the 
wind parameters, we use = 2.6 V esc and M predicted by the Vink et al. (2000) formula. 
The radially dependent wind density is determined from the mass conservation equation 
using a /3-velocity law V(r) = i>oo(l — -R*/r) /3 assuming a /3 = 0.8 (Pauldrach et al. 1986; 
Miiller & Vink 2008). The wind electron density is derived assuming hydrogen and helium 
are fully ionized. 

The energy dependent mean intensity (erg cm -2 s" 1 eV~ x str -1 ) used to calculate 
the photoionization rates has two dominant contributions: 1) the photospheric radiation 
field for a given T e ff and appropriate log g (using the TLUSTY grid of models from OS- 
TAR2002, Lanz & Hubeny 2003), along with the standard geometric dilution factor, and; 2) 
the XUV+X-ray radiation field, Jx- This Jx is specified by two parameters, the hot plasma 
temperature, Tx, and the column emission measure, AEMq (cm -5 ) such that Jx{E,Tx) = 
AEMc e(E,Tx) where e is the energy-dependent emissivity (erg cm 3 s" 1 eV _1 str -1 ) taken 
from the APED data. The basic assumption is that Jx represents an "effective" XUV+X-ray 
mean intensity at each given radial location, i.e., the wind contains a finite but small level 
of XUV+X-ray radiation distributed throughout the wind which seems to be supported by 
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observations as discussed in §2. All calculations presented in this section use a T x = 1 MK 
so we can examine the maximal effects of XUV radiation on q (P v) and g (S v). The total 
input mean intensity is determined for an energy grid from 8 eV to 2 keV. 

Since we are concerned with studying ionization effects for all luminosity classes over a 
large range in T e ff for which there is P V data, it is advantageous to define a new parameter 
that is dependent on both X-ray and stellar parameters. By defining Fx as the energy 
integral of 47rJx(E, T x ) above the He II edge (erg cm~ 2 s _1 ), then our fundamental adjustable 
parameter used in this study is defined as Fx / F* where F* is the total stellar photospheric 
flux (LboI/ '47T.R* 2 ). Hence, for a given F*, Fx/F*, and Tk, AEMq can be extracted directly 
from 



where A(Tk) (erg cm 3 s -1 ) is the total energy integral of 4-7re(T,Tx) above the He II edge. 
Note that Fx / F* is not the same as the well known observed X-ray to bolometric luminosity 
ratio, Lx/ L-Qoi, because Tx in this ratio is an "observed" quantity, i.e., a measure of only 
those X-rays capable of escaping the stellar wind, and our Tx is defined as an intrinsic total 
mean intensity (in flux units) where the majority of this emission (i.e., XUV) resides in an 
observational window that is inaccessible due to wind and ISM attenuation. 

The predicted q (P v) dependence on Fx / T* for supergiants, giants, and main sequence 
stars is shown in Figure 4. Also shown in this Figure are the data points from F06 that 
correspond to the required q (P v) values if all stars have their traditional M which we will 
use to determine the constraints on Fx/F*. This deficit in q (P v) relative to unity led F06 
to conclude that M needs to be reduced. Figure 4 shows a strong dependence of q (P v) 
on F x /F*, and indicates that a range in F x /F* between (0.3 — 10) x 10~ 7 can explain 
the observed q (P v) for all luminosity classes. For example, from Figure 4, the observed 
q (P v) of the four supergiants at T e ff = 35000 indicate a Fx/F* ps 2.5 x 10~ 7 . This implies 
a XUV+X-ray flux w 2 x 10 7 erg cm -2 s" 1 and a AEM C ~ 2.5 x 10 29 cm -5 (using Eq. 
1). From X-ray analyses of OB stars we cannot directly determine AEMc since only the 
volume emission measure EMy (cm -3 ) can be deduced from observations. If we assume that 
the XUV+X-ray radiation arises from a spherically shell at the assumed radius (r = 1.7 R* 
where R* = 20.5R Q ), then the "intrinsic" EMy is w 1.8 x 10 55 (using 4vrr 2 AEM C ) which 
is comparable to the lowest energy line "observed" EMy (ps 5 x 10 55 for N vil) derived from 
Chandra HETGS observations (e.g., Wojdowski & Schulz 2005). Since the intrinsic EMy is 
expected to be > the observed EMy (optical depth effects), the XUV+X-ray flux required 
to reduce q (P v) is well within the observational limits. 

Now we examine the effects of XUV+X-ray radiation on S v by considering the ratio 
q (P v)/q (S v). As shown in Figure 5, for the case when F x = 0, this ratio is ps 1 for a 
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large range in T e ff which supports the F06 S V surrogate argument. However, as F x jF* 
increases, g(P v)/g (S v) decreases which means that q (S v) is significantly less sensitive 
to the XUV+X-ray radiation as compared to the dependence of q (P v). This is a direct 
consequence of the difference in the P v and S V photoionization rates shown in Figure 3 
(see §2), and invalidates the sulfur surrogate argument. In general, for all luminosity classes, 
g(P v)/g (S v) reaches a minimum value of ~ 0.14 over most of the O-star spectral range. 
The required Fx/F* to produce this minimum value is dependent on luminosity class as 
shown in Figure 5. These results imply that the M derived from FUSE observations could 
be underestimated by almost an order of magnitude. 



4. Discussion 

We have demonstrated that the presence of XUV radiation embedded throughout a 
stellar wind in sufficient amounts (well within the observational constraints) can lead to a 
significant depletion of q (P v) by using only outer shell photoionization processes. There- 
fore, the discrepancy between P v derived M with those obtained from density-squared 
diagnostics is far less severe than suggested by F06. Hence, these stars can have M that are 
again back within a factor of 2 range of the traditional M (consistent with clumped wind 
predictions). This also means that the XODUR does not require an alternative explanation. 

Gudel & Naze (2009) mentioned that there may be a problem with XODUR since the 
wind opacities used by Waldron & Cassinelli (2007) were determined from a more highly 
ionized wind. Although it is true that the opacity at low energies can be sensitive to the 
assumed wind ionization structure, the key point is that the XODUR is based entirely on 
emission lines at energies > 0.56 keV where the opacity is almost identical to the "cold" 
ISM opacity above this energy (e.g., see Fig. 2 of Waldron et al. 1998), and the ISM opacity 
represents the unsurpassable upper limit to any wind opacity. Therefore, regardless of the 
wind ionization state, the XODUR is certainly valid for those radii derived from the He-like 
fir lines of Ne IX, Mg XI, Si xill, and S XV, and radii derived from O VII may be marginally 
dependent on the wind ionization state. 

We have shown that the inclusion of XUV radiation can resolve the M discrepancy, 
and there is observational support that this XUV radiation must be distributed throughout 
these winds based on the conclusion of M03 regarding C V. But, what is the source of this 
emission? Since these winds are clumped, this emission is likely to originate from bow shocks 
forming around clumps. The Cassinelli et al. (2008) wind bow shock model predicts that the 
temperature dependence of the emission measure scales as (T/T max )~ 4 ^ 3 , where T max is the 
temperature at the apex of the bow shock. Hence, there should be a considerable amount 
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of XUV radiation, even for rather strong shocks, located at any radius. Although this XUV 
radiation is unobservable since essentially all of this radiation will be absorbed by the stellar 
wind and ISM, the strength of this emission can be determined indirectly from analyses of 
lower energy spectral bands as demonstrated in this paper. 

Future studies now need to focus on alternate explanations as to why the X-ray lines are 
symmetric and nearly un-shifted. Several possibilities have been proposed: wind porosity 
effects (Oskinova et al. 2004, 2006); asymmetric mass outflows (Mullan & Waldron 2006), 
and; bow shocks forming around wind clumps (Cassinelli et al. 2008) and stellar ejected 
plasmoids (Waldron & Cassinelli 2009). In addition, Oskinova et al. (2007) claim that the 
X-ray line symmetry and P V problems can both be resolved, with no need to reduce M, using 
their planar shock model, porosity effects, and macro-clumping. However, we argue that bow 
shocks will undoubtedly form around these macro-clumps and the impact of the resultant 
excess XUV radiation on q (P v) needs to be examined. Regardless as to which explanation 
is applicable, our work has re-emphasized the importance of including the full XUV+X-ray 
energy spectrum when exploring the effects of radiation on stellar wind ionization structures. 

This work was supported by NASA ATFP award NNH09CF39C. 
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Fig. 1. — Comparison of radiative cooling curves determined from the RS, MEKAL, and 
APED emissivites. The dashed-line curves represent the contribution of the XUV radiation 
(54.4 - 124 eV) to the total cooling curve. The solid blue line represents the contribution 
from the important collection of XUV emission lines between 65 and 73 eV (See §2 and Fig. 
2). 
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Fig. 2. — Comparison of the energy dependent X-ray emissivity (Tx = 1 MK) with the 
photoionization energy ranges of the relevant ions of carbon, phosphorous, and sulfur. This 
shows the strong collection of XUV emission lines just below the S V threshold energy at 
72.68 eV, and above this energy, the XUV radiation is essentially devoid of emission lines 
(see §2). The four dominant XUV lines are Fe IX (72.47 eV), Fe X (71.04, 69.95 eV), and 
Fe XI (68.72 eV). The vertical dashed line is the He II edge (54.4 eV). 
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Fig. 3. — The XUV+X-ray photoionization rates for C IV, P v, and S v as a function of 
temperature assuming that 4tt times the total energy integrated XUV+X-ray mean intensity 
above 54.4 eV is fixed at 10 8 erg cm -2 s _1 for each temperature. The dashed-line curves 
show these rates as determined by using only the X-ray energy band (> 124 eV). 
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Fig. 4. — The effects of XUV radiation on q (P v) as a function of T e ff for different values 
of log (Fx/F*) for the three luminosity classes. The black dashed-line represents the Fx = 
case. The filled circles are the observationally determined q (P v) assuming traditional M 
from F06 (see §3). 



- 15 - 




30000 35000 40000 45000 

T e ff 

Fig. 5. — The ratio of q (P v)/q (S v) as a function of T e f f for different values of log (Fx / -F*) 
for the three luminosity classes. The black dashed-line represents the Fx = case and this 
ratio is ~ 1 for T e ff < 39000 K. See discussion in Section 3. 



